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LIKELY VALUES OF THE COSMOLOGICAL CONSTANT 

Hugo MarteP’^, Paul R. Shapiro^’^, and Steven Weinberg‘S’® 

ABSTRACT 

In theories in which the cosmological constant takes a variety of values in 
different “subuniverses,” the probability distribution of its observed values is 
conditioned by the requirement that there must be someone to measure it. This 
probability is proportional to the fraction of matter which is destined to condense 
out of the background into mass concentrations large enough to form observers. 

We calculate this “collapsed fraction” by a simple, pressure-free, spherically sym¬ 
metric, nonlinear model for the growth of density fluctuations in a flat universe 
with arbitrary value of the cosmological constant, applied in a statistical way to 
the observed spectrum of density fluctuations at recombination. From this, the 
probability distribution for the vacuum energy density pv for Gaussian random 
density fluctuations is derived analytically. (The conventional quantity Aq is the 
vacuum energy density in units of the critical density at present, Aq = pvjPent, 
where Pcrit,o = ^Hq ISttG.) It is shown that the results depend on only one quan¬ 
tity, cr^p, where cr^ and p are the variance and mean value of the fluctuating 
matter density held at recombination, respectively. To calculate cr, we adopt 
the hat CDM model with nonzero cosmological constant and hx the amplitude 
and shape of the primordial power spectrum in accordance with data on cosmic 
microwave background anisotropy from the CODE satellite DMR experiment. A 
comparison of the results of this calculation of the likely values of pv with present 
observational bounds on the cosmological constant indicates that the small, pos¬ 
itive value of Pv (up to 3 times greater than the present cosmic mass density) 
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suggested recently by several lines of evidence is a reasonably likely value to 
observe, even if all values of pv are equally likely a priori. 


Subject headings: cosmology: theory — galaxies: formation 


1. INTRODUCTION 

Though the evidence is still equivocal, there are persistent hints that the vacuum energy 
density®/?!/ is positive, and up to 3 times greater than the present cosmic mass density 
From the point of view of fundamental physics, such a value seems absurd. Crude estimates 
indicate a value of pv some 120 orders of magnitude greater than po, and while it is hard 
enough to imagine any sort of symmetry or adjustment mechanism that could make pv 
vanish (for a litany of failed attempts, see Weinberg [1989]), it would be even more peculiar 
for fundamental physical theory to dictate a non-zero value for pv that happens to be 
comparable to the cosmic mass density po at this particular moment in the history of the 
universe. 

As far as we know, the only way to understand a value of pv comparable to po is based on 
a weak form of the anthropic principle. In several current theories the cosmological constant 
does not have a hxed value, but takes a variety of values with varying probabilities. For 
instance. Hawking (1983, 1984) showed that the introduction of a three-form gauge held 
yields a state vector for the universe that is a superposition of terms with different 
values for the cosmological constant. Coleman (1988a) subsequently showed that the effect 
of wormholes in quantum gravity is to make the state vector a superposition of terms, in 
which any coupling coefficient in the Lagrangian that is not hxed by symmetries takes all 
possible values.® Also, in chaotic inhation (Linde 1986, 1987, 1988) the observed big bang is 


®By Pv is meant the sum of a term A/SttG, where A is the cosmological constant appearing in the 
Einstein held equations, plus the contribution to the vacuum energy density from quantum huctuations.The 
conventional quantity Aq is the vacuum energy density in units of the critical density at present, Aq = 
Pv/Pcrit,o, where pcnt.o = 3H^/SttG. 

^For a review and earlier references, see Ostriker & Steinhardt (1995). This conclusion has been recently 
challenged by preliminary results of measurements by Perlmutter et al. (1996) of redshifts and distances for 
distant Type la supernovae. 

®Coleman (1988b) subsequently concluded that the probability distribution of pv is sharply peaked at 
zero, as had previously been argued by Hawking (1983, 1984) and Baum (1984). This conclusion has been 
challenged by Fischler et al. (1989), and it will be assumed here that there is no sharp peak at py = 0. 
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just one of an infinite number of expanding regions, in each of which the various holds that 
affect the vacuum energy can take different values. For brevity we will refer to parts of the 
“universe” in which the cosmological constant takes different values, such as terms in the 
state vector, local bangs, or whatever, as subuniverses. 

In any theory of this general sort the measured value of the vacuum energy density 
pv would be much smaller than the value expected on dimensional grounds in elementary 
particle physics, not because there is any physical principle that makes it small in all sub¬ 
universes, but because it is only in the subuniverses where it is small that there would be 
anyone to measure it. This paper will show how to calculate the probability distribution of 
the values of pv that would be observed under these circumstances. 

An earlier paper (Weinberg 1987) pointed out that the anthropic limit on the value of 
Pv for py > 0 arises from the requirement that pv should not be so large as to prevent 
the formation of galaxies. This paper suggested that this requirement implies a value of 
Pv roughly comparable to the cosmic density of nonrelativistic matter at the time that the 
earliest galaxies form, because, if pv were much larger than this, then galaxies could not 
form and there would be no observers, while there did not seem to be any reason for pv 
to be much smaller than this. Since then galaxies have continued to be found at higher 
and higher redshifts, and hence at higher and higher values of the cosmic mass density, and 
it is becoming clear that such values of pv are already ruled out. A galaxy with redshift 
z ^ 4 was formed when the cosmic mass density was more than (and perhaps considerably 
more than) (1 + z)^ sa 125 times the present mass density, which is much greater than the 
observational upper limit on pv- 

On the other hand, it is much more likely that the value of pv in our subuniverse is 
comparable to the average or median value measured by astronomers in all subuniverses, 
rather than the anthropic upper bound, so that its value should be compared with the 
cosmic mass density at the time of formation of typical galaxies, rather than of the earliest 
galaxies.® Here we will present a detailed Bayesian analysis, which allows a calculation of 
the probability distribution of pv from a knowledge of the spectrum of density fluctuations 
at recombination. The results suggest a much smaller likely value of pv than the anthropic 


®See Weinberg (1996). This is essentially the same as what was called the “principle of mediocrity” by 
Vilenkin (1995a, 1995b, 1995c, 1995d). Vilenkin did not undertake a detailed calculation of the probability 
distribution of the cosmological constant. A calculation of this sort was done by Efstathiou (1995), but 
it contained some errors (Vilenkin 1995d; Weinberg 1996). Efstathiou’s calculation was done numerically, 
using linear perturbation theory and what is believed to be a realistic model of initial perturbations, while 
the calculation presented here is thoroughly nonlinear, but concentrates on a single spherically symmetric 
density fluctuation, so that it is possible to understand the results analytically. 
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upper bound, a value that may not be in conflict with present observational bounds. 

In §2 we describe how to calculate the probability distribution for the vacuum energy 
density pv that would be observed in various subuniverses in which there is someone to 
observe it. For pv > 0, this is simply related to the fraction of matter that condenses into 
galaxies. We evaluate this fraction and the resulting probability distribution in §3 under the 
assumption of Gaussian random density fluctuations in the cosmic mass density at the time 
of recombination. The results depend only on the standard deviation a of these fluctuations 
at the time of recombination. In §4 we calculate cr, adopting the cold dark matter model for 
the power spectrum of these density fluctuations, and assuming a flat universe with nonzero 
cosmological constant (sometimes referred to as the flat “ACDM” model). The amplitude 
of the density fluctuations is hxed by the data on cosmic microwave background anisotropy 
from the CODE satellite DMR experiment. The results of this calculation of the likely values 
of Pv are presented in §5, and compared with the range of values of the cosmological constant 
allowed by current observational and theoretical constraints. A summary and conclusions 
are presented in §6. 


2. THE PROBABILITY DISTRIBUTION 


We assume that the a priori probability of a net vacuum energy density between pv 
and Pv + dpv is 'P{pv) dpv , where V{pv) is some smooth function of pv, with no special 
behavior near pv = 0. What we want is the probability distribution 'Pohs{pv) that a random 
observer in any subuniverse will measure values of pv in a given range. According to the 
principles of Bayesian statistics, this is given by 


'Pohs{pv) 


/o“ A{pv)V{pv) dpv 


( 1 ) 


where A{pv) is the mean number of astronomers making independent measurements of the 
vacuum energy density in subuniverses with vacuum energy density pvd^ 

In calculating the quantity (1), we note that the range of anthropically allowed values 
of Pv is so much smaller than the energy densities typical of elementary particle physics, 
that, within this narrow range, we can take the a priori probability distribution V{pv) to 


have not thought through the problems associated with inhnite subuniverses, where A if nonzero is 
inhnite. Presumably in this case we should take A to be the number of independent astronomers per same 
hxed number of baryons. 
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be constant. The value of this constant then cancels in equation (1), which becomes 

A{pv) 


Vohsipv) = 


fcT A(pv)dpv 


( 2 ) 


The evolution of galaxies and astronomers depends on a variety of constants of nature 
other than py, and the values of these other constants in the various subuniverses may be 
correlated with the values of py, but the range of values of pv that are anthropically allowed 
is so small compared with the energy densities typical of elementary particle physics that, 
within this range, we can take all other fundamental constants to have hxed values, the values 
we observe in our subuniverse.Also, once a fluctuation in the cosmic mass distribution 
undergoes gravitational condensation, its subsequent evolution is essentially independent of 
py, so the ratio of astronomers to mass in galaxies may be taken as independent of py. The 
number of astronomers ^(py) who can measure py in any subuniverse should therefore be 
proportional to the fraction jF(py) of matter incorporated in galaxies, so that equation (2) 
may be written 

Hpv) 


Vohsipv) = 


/o“ :f{pv) dpv 


( 3 ) 


To calculate jF(py), we note that the spectrum of initial density fluctuations at re¬ 
combination can be regarded as independent of py, because values of py for which galaxy 
formation is possible are much smaller than the cosmic mass density at or before the time 
of recombination. Similarly, it is reasonable to suppose that the total amount of matter in 
a subuniverse in theories of chaotic inflation is independent of py within the narrow range 
of values of py that are anthropically allowed. Our problem is then to calculate the fraction 


might be asked why within this range the a prion distribution V{pv) is not a power of py rather 
than a constant? A power law would mean that there is something special about the value py = 0, since a 
power law distribution function would have to vanish or blow up there. But the essence of the cosmological 
constant problem is that we do not know of anything in fundamental physics that gives a special signihcance 
to the value py = 0, which requires a precise cancellation of a coupling coefhcient in the Lagrangian by 
radiative corrections. (Analogously, although the probability distribution of temperatures in the Antarctic 
ice must vanish at zero degrees Kelvin and zero degrees Centigrade, since these are the limits of the range 
of temperatures in which water freezes, we would hardly expect it to vanish or diverge at zero degrees 
Fahrenheit.) The same argument applies to a power-law dependence on logpy. Of course we are not 
assuming that V{pv) does not contain terms that vary as positive powers of py, but only that there is also 
a constant term, which then naturally dominates for the very small values of py that are consistent with the 
appearance of astronomers who can measure py. 

^^If various constants of nature and initial conditions vary from one subuniverse to another independently 
of the values of py, then Vohs{pv)dpv is the probability that, */the other constants and initial conditions 
take the values we observe, then the vacuum energy density will be observed to be between py and py -\-dpv ■ 
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T of matter that undergoes gravitational condensation into galaxies as a function of /jy, 
for fixed initial conditions at recombination. Ironically, while the tininess of observationally 
allowed values of pv creates the cosmological constant problem in the hrst place, it is the 
tininess of the range of anthropically allowed values of pv that offers the possibility of a 
realistic calculation of Vohs{pv)- 

To see how this can work in practice, we will carry out an illustrative calculation of J^[pv) 
and use it to calculate both the integrated probability distribution and the mean and median 
values of pv observed by all astronomers in the subuniverses that contain astronomers. 
Earlier work (Weinberg 1987) used a very simple model (Peebles 1967) of galaxy formation 
from isolated spherically symmetric pressureless fluctuations. This calculation was improved 
in the report of a recent conference talk (Weinberg 1996), by using the well-known model of 
Gunn & Gott (1972), which also assumes isolated fluctuations with spherical symmetry and 
zero pressure, but includes the infall of matter from outside the initially overdense ball. Here 
we will also take into account the facts that, with space hlled with fluctuations, there is a 
limit to the mass that can accrete onto any one fluctuation, and that there must be regions 
of negative as well as positive overdensity. 


Consider a spherically symmetric pressureless fluctuation, consisting at recombination of 
a spherical core of volume V and positive average fractional overdensity 6 [i.e. 6 = [p — p)lp^ 
where p is the average density inside V and p is the cosmic mean density at recombination], 
surrounded by a spherical shell of volume U of constant fractional underdensity, taken to 
have the value [VlU)6^ so that the average overdensity within the whole fluctuation is zero. 
Outside this shell are other fluctuations, about which we do not need to say anything, except 
to assume that they do not seriously interfere with the spherical symmetry of the fluctuation 
in question. For simplicity, we will take VjU to be the same for all fluctuations. 


ifl3 


As shown in earlier work (Weinberg 1987) the core will undergo gravitational collapse 


8 > 


mpvV''' 

500p ) 


(4) 


In addition, that portion U' of the outer shell will fall into the core, for which the average 
fractional overdensity within the volume V -\- U' saturates this inequality: 


V8 - U'{VIU)S _ / 729pv y^^ 
V + U' \ 500p ) 


(5) 


was originally assumed (Weinberg 1987) that the overdensity 6p was uniform, but these results actually 
hold for arbitrary spherically symmetric fluctuations, with 6 interpreted as the average fractional overdensity. 



The fraction of the total mass piU + V) that suffers gravitational contraction will be 

(1 + h) y + [1 - {viu) 8 ] u' 


HS.py) = 


W + V) 


( 6 ) 


Solving equation (5) for U\ we hud for the fraction of mass that undergoes gravitational 
contraction 

1 + (729/jv/500^)^/3 


[{729pv/500py/^ + {V/U)6\ 
We require that the total density be everywhere nonnegative, so that 

S < (U/V) . 


(7) 


( 8 ) 


For 6 satisfying this inequality, equation (7) gives iF < 1, so that no fluctuation can get more 
than its fair share of mass. 

In what follows we will assume that the fluctuation number density A/'(h) is negligible 
for initial fluctuations that are not everywhere weak, so that we will be integrating only over 
fluctuations for which h ^ 1 and 6 U/V. Also, for any anthropically allowed cosmological 
constant, pv is much less than the mass density p at recombination, so we will drop the 
term (729pv/500p)^/^ in the numerator (but not the denominator) of the fraction in square 
brackets in equation (7). The fraction of mass winding up in galaxies is then 


Hpy) = I 


(729py/500p)l/3 


dS A'(h) .F(h, pv) 


dS 


{VIU)8M{6) 


(729py/5oop)i/3 (729pv/500p)^/^ + (y/f/)h ’ 


(9) 


where J\f (8) d8 is the fraction of all positive fluctuations that have average core fractional 
overdensity between 8 and 8 + d8, normalized so that 

roo 

/ A'(h) d8=l . 


The normalization integral in equation (3) can be calculated by interchanging the order 
of integration over 8 and py, and expressing pv in terms of a dimensionless variable x dehned 
by 


pv = 


729 


Equation (9) then gives 


roo '^00 

/ Hpv) dpv = (V) (C/fO/o(vyto , 


( 10 ) 



- 8 - 


where Io{s) is the function 


dx 1 ,, 

Io(s) = / - =- s-s^ln, 

Jo X s 2 Vl + s 


and the brackets denote an average over all positive fluctuations 


{/(«))=/ m)f(s)ds. 

Jo 


( 11 ) 


( 12 ) 


The normalized probability distribution (3) for the observed vacuum energy density is then 


Vohsipv) 


243 1 

500 {63)Io{VlU) 

r _ 

d(729py/500p)l/3 (729pv/500p)^/^ + (y/f/)h ’ 


(13) 


with all quantities on the right-hand side referring to the time of recombination. 


In using equation (13) we will need to make some assumption about the shape parameter 
s = V/U. The value s = 0 corresponds to the limit in which each positive fluctuation is 
isolated, surrounded by an inhnite volume of compensating underdensity (at a total density 
arbitrarily close to its mean value p), the case considered by Weinberg (1996). Values of s 
much greater than unity correspond to the limit in which the additional mass associated 
with the compensating underdense volume U is insignihcant compared with that contained 
within the positive fluctuation in volume V, the case considered in Weinberg (1987). The 
value s = 1 corresponds to the case in which every positive fluctuation is surrounded by 
an equal volume of compensating negative fluctuation. This latter value is the one most 
relevant to a Gaussian-random distribution of linear density fluctuations, since the volumes 
occupied by positive and negative density fluctuations of equal amplitude are exactly equal 
in that case. Thus we will concentrate on the value s = 1 when we apply our analysis to the 
observed universe in what follows. Fortunately, as we shall see, most of our results will turn 
out to be almost independent of the value chosen for s. 


Strictly speaking, equation (13) times dpv gives the probability that, if the vacuum 
energy density is positive, then it will be observed to be between pv and pv + dpv- For 
pv < 0, the anthropic bound on pv is set by the condition (Barrow & Tipler 1986) that 
the subuniverse should survive long enough for intelligent life to arise. It is plausible that 
Pv > 0 is strongly favored anthropically (Weinberg 1996) as well as observationally, but we 
will make no attempt here to calculate the probability distribution for negative values of pv- 


Fquation (13) can be used in various ways. One is to calculate the mean of various 
powers of pv - Carrying out the same exchange of integrations and change of variables as in 
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our calculation of equation (10), we find 


<' >= V‘ LAViU) 

V729/ («3) laiVIU) 


(14) 


where L is the function 


n dx 


Us) = / — ■ ( 15 ) 

t/o X ~Y S 

The average <Py> in equation (14) is taken over all subuniverses, not over fluctuations as 
in (h), and again all quantities on the right-hand side are to be evaluated at recombination. 
In particular, the mean value observed for pv (if pv is positive) is 


with 


500 {6^)pIi{V/U) 

m (h3) loiv/u) ’ 


^i(^) 


1 

5 



+ 3“^ + s® In 


s 

1 + s 


(16) 


( 17 ) 


Fortunately the ratio IiiVjU) jIo{VjU) in equation (16) turns out to be nearly constant; 
it drops from a value 0.5 when s = VjU 1, corresponding to no infall, to a value 0.4 when 
s ^ 1, corresponding to well separated fluctuations. It is therefore not very important what 
value we choose for s. As mentioned earlier, it seems reasonable to assume the intermediate 
case s = 1, where overdense and underdense regions have equal volume at recombination; in 
this case, the ratio /i//o takes the value 

^ - ln2 

/i(l)//o(l) = = 0.46693 . 


The insensitivity of our results to the value of s suggests that they also may not be much 
affected by the crudeness of our treatment of the effect of one fluctuation on another. 

We could also use equation (13) to calculate the integrated probability 'P(> pv) that 
the vacuum energy density is greater than pv'- 

roo 

V{>pv)= V.Up'v)dp'v (18) 

■j Pv 


With the same reversal of integrations over p'y and 8 and the same change of variables as 
before, we hnd 


V{> pv) 


(h3 /((729pv/500p)i/3/h , V/U) ) 

W)WJu) 


(19) 
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where for t < 1, 

j -^^dx = hl-f)-s{l-t) + s‘^ln(^-^'], ( 20 ) 

Jt X s Z \t s / 

and /(f,s) = 0 for f > 1. As in our calculation of < pv >, this is insensitive to the precise 
value of s = V/U. For very small py, equation (19) of course approaches unity, whatever 
the value of s. For very large py, the only fluctuations that contribute within the integral 
over 6 in equation (19) are those with 6 near the lower limit of the integral, for which the 
lower limit of the integral (20) is near unity, where this integral behaves as 1/(1 + s). But 
then V{> Pv) has an s-dependence proportional to [/o(s)(l + s)]“^, which only rises from 2 
to 3 as s rises from zero to inhnity. 


3. Gaussian Density Fluctuations 


To go further, we must make some assumption about the form of the fluctuation prob¬ 
ability distribution ^”((5) at an early epoch, such as that of recombination. Current data 
on the anisotropy of the cosmic microwave background on large angular scales and on the 
large-scale clustering properties of galaxies, as well as theoretical predictions of the origin 
of density fluctuations by quantum processes in the early universe in inflationary cosmology 
models, are consistent with the assumption that the primordial fluctuations were isotropic, 
Gaussian random noise of very small amplitude. In this case the fluctuation distribution has 
the form 



( 21 ) 


The mean values of powers of 6 are given in terms of the variance X by 

2JV/2^JV m 

) ■ 


( 22 ) 


Fquation (14) then gives 

<i>y>=r(; 


3n + 4^ /lOOO • 2i/VVy 4(s) 

V 729 ) 7o(iy’ 


and in particular, the mean vacuum energy density is 


625(27r)^/^cr®p 



(23) 


< pv > 


243 


(24) 



where, as before, s = V/U. This gives the numerical values 


< Pv >= X < 


f 2.5788 s = 0 , 

3.0103 s = 1 , 

I 3.2235 s = oo 


( 25 ) 


Also, using the Gaussian distribution (21) in equation (19) and writing 

~ V 5ooh ) U; ’ 


we hud the differential probability distribution 

e~^ dx 


'Pohs{pv) dpv = 


2Io{s) Jfj -|- ^ 1/2 ’ 


where 


/3 = 


1 


|^ 729pv \ 


2/3 


2(72 500p ) 

The probability of a vacuum energy density greater than pv is then 

1 


iP(>pv) = (l+/3)e-^- 


2/o(s) 


e In 


(26) 


(27) 


s \x/ 


>dx . (28) 


By combining equations (24) and (27), we see that the parameter fi in equations (26) and 
(28) may be expressed in terms of s = VjU and the ratio pvl < Pv> of the vacuum energy 
density to its mean value 


/3 



Pv h{s) 
< Pv > Io{s) 


2/3 


(29) 


Thus the probability of observing a vacuum energy density in a certain range depends only 
on the values of pv I < Pv> and s. The analysis of initial fluctuations will enter here only as 
a means of calculating the parameter pcr^, and hence <pv>- 

In using these formulas, it is important to note that during the era of recombination, 
when Pv is negligible and fluctuations are small, a grows as /2/3 and p falls as ^ so the com¬ 
bination (T^p is time independent. Therefore equations (24)-(28) show that < py >, 'Pobs(/0v), 
and P{> Pv) do not depend on what we take as the precise moment of recombination at 
which a and p are evaluated. 
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We have plotted the diflerential probability dVjda = < pv > 'Pohs{pv) versus 

a = Pv I < Pv > in Figure 1, for a range of values of s. This hgure shows that the probabil¬ 
ity distribution drops exponentially for large py, and is remarkably insensitive to the value 
of s for all py, except for py ^ 1. We have also plotted the differential probability per 
logarithmic interval of a, a'Pobs(<n), versus a [i.e. ddPjdAoga = pv'Pohs{pv)\ in the bottom 
panel of Figure 1. This quantity is peaked at values of a = 0.7 — 0.8, almost independent 
of s. Figure 2 shows the integrated probability 'P(> py) versus the dimensionless quantity /3 
dehned by equation (27), for various values of s. Again, we see that this probability is only 
weakly dependent on s. We can evaluate V{> pv) exactly in a few special cases 


V{> pv) 




Ip-P 


s = 0, 

s = oo & /3 1 . 


(30) 


Fquation (28) may be used to calculate the median value (py)i/ 2 , for which 
'P[> (py)i/ 2 ] = 1/2. For isolated fluctuations s = 0, and the result is 


(dn)i /2 = 


4(ln2)3A 


<PV>= 0.43411 <py>= 1.1195 (jV 


(31) 


For other values of s the median must be calculated numerically, with the result that 

' 1.1195 
1.4851 


{Pv)l/2 = <pv> X \ 


0.43411 

0.49335 = a^p 

0.51480 


s = 0 , 
s = 1, 


(32) 


1.6595 


s = oo 


Again we see the insensitivity of our results to the shape parameter s = V/U. 


It is also interesting to ask what is the range of reasonably likely values of py? For 
instance, what is the range of values of pv for which only 5% of astronomers in all subuni¬ 
verses would observe smaller values, and only 5% would observe larger values? By setting 
'P(> py) in equation (28) equal to 0.95 and 0.05 and solving for py, we hnd that this range 
has the lower bound 


(/5v)o.95 


< py > X < 


0.00874 

0.01959 

0.02359 


(J^p X 


0.02254 

0.05897 

0.07604 


s = 0, 

■5 = 1, 

S = (DO , 


and the upper bound 


= <py> X < 


3.9005 

3.6914 

3.6157 


a^p X 


10.0586 

11.1122 

11.6552 


s = 0 , 
■5 = 1, 


(33) 


(/5v)o.05 


OO . 


(34) 











Fig. 2.— Integrated probability 'P(> fi) versus /3, for s = 0.01, 0.1, 1, 10, and 100. 

Although the lower bound (33) is evidently somewhat sensitive to the shape parameter s, we 
see that for all values of s the distribution of pv values is quite broad; it would not be very 
unlikely for a subuniverse to have a value of pv that is 50 times smaller or 3.7 times larger 
than the average. On the other hand, it would be extremely unlikely to observe a value of 
Pv which differs from the mean by more than a few orders of magnitude. Not only are large 
values of pv unlikely, therefore, as we might previously have guessed based upon the fact 
that large pv suppresses galaxy formation, but values of pv extremely close to zero are also 
unlikely; there are simply too many other subuniverses to observe which have larger values 
of Pv but not large enough to prevent galaxy formation. 

For some purposes it is useful to have an analytic ht to the integrated probability. The 




following generally works quite well: 


(3s + 2)e ^ — se 
2s + 2 


(35) 


VI, > pv ) 


Notice that this ht is consistent with the special cases mentioned above. The absolute 
difference [V — f\ and the relative difference liV — f)/V\ are always less than 0.025 and 0.028, 
respectively, for 0 < s < 1, while these errors increase to 0.027 and 0.092, respectively, for 
1 < s < (DO. Because the approximation formula (35) gives V{> pv) as a quadratic function 
of e“^, it is easy to solve equation (35) for the value of [3 and hence of pv/ < pv> that give 
any specihc value for 'P(> pv)- For instance, for s = 1 equation (35) yields a median vacuum 
density {pv)il 2 = 0.482 < pv >, as compared with an exact value {pv)il 2 = 0.493 < pv >, 
while for s = oo equation (35) yields {pv)i /2 = 0.568 < pv >, as compared with an exact 
value (/ 0 v)i /2 = 0.515 <pv>- For s = 0, equation (35) is exact. 

Finally, it is also useful to compute the probability that a is observed to be within the 
interval [q(i,Q( 2 ], according to V{ai < a < 0 : 2 ) = V{> cti) — V{> 0 : 2 ). We present these 
interval probability results in Table 1, for s = 1. We also plot in Figure 3 the interval 
probability isocontours such that, for any given value of a, call it a*, the curves show the 
values of a above and below this a* at which the interval probability has some particular 
value, as labelled. In particular, each curve in the (a*, Q()-plane in Figure 3 corresponds to 
the locus of points which satisfy the equation V{> a*) — V{> a) = constant, if a > a*, 
or 'P(> a*) — V{> a) = constant, if a < a*, where the labels indicate the values of the 
constant. 


4. Evaluation of a 

4.1. Filtered Density Fluctuation Spectrum 

The mean value < py > as well as the probability distribution 'Pohs{pv) and the in¬ 
tegrated probability 'P(> pv) have been expressed in equations (24)- (29) in terms of the 
variance in the fluctuation distribution (21). Now we must consider how to calculate a. 

From equation (22), we have = (h^). But the variance which is appropriate for our 
purpose here, is that which reflects the range of wavenumbers which might possibly contribute 
to the formation of gravitational condensations that are large enough to lead to “astronomer 
formation.” Only wavenumbers corresponding to density fluctuations encompassing such 
sufficiently large masses should be allowed to contribute. This implies that the appropriate 
a for our purpose here is one calculated by hltering the underlying density held to eliminate 
the contribution from small wavelengths. This is accomplished by smoothing the density 
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Table 1. THE PROBABILITY THAT pv/ < /Jy > IS BETWEEN TWO VALUES^ 


ai 

0(2 

V{> ai) - V{> a2) 

ai 

0(2 

V{> ai) - v{> a 2 ) 

0.000 

0.001 

0.00355 

1.5 

1.6 

0.01499 

0.001 

0.002 

0.00317 

1.6 

1.7 

0.01366 

0.002 

0.003 

0.00299 

1.7 

1.8 

0.01245 

0.003 

0.004 

0.00287 

1.8 

1.9 

0.01140 

0.004 

0.005 

0.00277 

1.9 

2.0 

0.01044 

0.005 

0.006 

0.00269 

2.0 

2.1 

0.00959 

0.006 

0.007 

0.00262 

2.1 

2.2 

0.00882 

0.007 

0.008 

0.00256 

2.2 

2.3 

0.00813 

0.008 

0.009 

0.00251 

2.3 

2.4 

0.00750 

0.009 

0.010 

0.00247 

2.4 

2.5 

0.00693 




2.5 

2.6 

0.00642 

0.00 

0.01 

0.02820 

2.6 

2.7 

0.00585 

0.01 

0.02 

0.02268 

2.7 

2.8 

0.00561 

0.02 

0.03 

0.02023 

2.8 

2.9 

0.00512 

0.03 

0.04 

0.01858 

2.9 

3.0 

0.00477 

0.04 

0.05 

0.01731 




0.05 

0.06 

0.01629 

0.0 

0.5 

0.50350 

0.06 

0.07 

0.01543 

0.5 

1.0 

0.18691 

0.07 

0.08 

0.01468 

1.0 

1.5 

0.10271 

0.08 

0.09 

0.01404 

1.5 

2.0 

0.06294 

0.09 

0.10 

0.01346 

2.0 

2.5 

0.04097 




2.5 

3.0 

0.02777 

0.0 

0.1 

0.18090 

3.0 

3.5 

0.01937 

0.1 

0.2 

0.11175 

3.5 

4.0 

0.01382 

0.2 

0.3 

0.08510 

4.0 

4.5 

0.01004 

0.3 

0.4 

0.06863 

4.5 

5.0 

0.00740 

0.4 

0.5 

0.05712 

5.0 

5.5 

0.00553 

0.5 

0.6 

0.04850 

5.5 

6.0 

0.00418 

0.6 

0.7 

0.04178 

6.0 

6.5 

0.00318 

0.7 

0.8 

0.03638 

6.5 

7.0 

0.00245 

0.8 

0.9 

0.03197 

7.0 

7.5 

0.00189 

0.9 

1.0 

0.02828 

7.5 

8.0 

0.00148 

1.0 

1.1 

0.02518 

8.0 

8.5 

0.00115 

1.1 

1.2 

0.02253 

8.5 

9.0 

0.00092 

1.2 

1.3 

0.02024 

9.0 

9.5 

0.00072 

1.3 

1.4 

0.01825 

9.5 

10.0 

0.00058 

1.4 

1.5 

0.01652 

10.0 

CXD 

0.00249 


^For s = 1, where a = pv /< pv > 
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Fig. 3.— Interval Probability Isocontours. Each curve in the (a*,a)-plane is the locus 
of points with a constant probability that a subuniverse is observed to have a value of 
py/ < py > in the range between a* and a (i.e. \V{> a*) — V{> a)| = constant, where 
each curve is labelled with the value of this constant.) 

held before calculating the variance cr^, according to 

a^ = (h^(r)), (36) 

where 

h(r) = y h(x)lF(x — r)d^a;. (37) 

Here IF is a smoothing “window function,” and x and r are co-moving coordinates, which 
following convention we shall normalize to give the proper distance at present. This yields 
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the following familiar expression for the variance cr^ 


a 


2 


1 


0 


P{h)W'^{kR)k'^dk, 


(38) 


where P{k) is the power spectrum (assuming statistical translation and rotation invariance) 

P(|k|) = J d^x {6{x + r)6{r)) (39) 

and W{kR) is the Fourier transform of the window function 

WikK) = I d\ , (40) 

with R a length parameter to be specihed below, introduced to make the argument of W 
dimensionless. (The bracket in eq. [39] implies an average over space; the assumptions of 
isotropy and homogeneity ensure that P depends only on k = |k|.) The window functions 
in which we are interested here are those which hlter out modes of wavelength smaller than 
some characteristic value R. There are two conventional choices for the window function, 
the Gaussian window function, 

WG{u) = e--"l\ (41) 


and the Top-Hat window function, 

- 3 

= —(sin u — u cos u) , 

7/3 


(42) 


(Peebles 1980). The baryonic mass associated with a density fluctuation of wavelength close 
to the hlter scale is given by 


M} = 


PbqR^ (Gaussian), 
(47r/3)pBo-R^ (Top-Hat), 


(43) 


with both R and the cosmic mean baryon density pBo evaluated at the present. The radii 
for which both window functions enclose the same mass are thus related by 


Rg 

i?TH 


(47r/3)^/3 

{2ttYP 


0.6431 


(44) 


We shall occasionally refer to a maximum wavenumber fc ma y 
ing to a wavelength i?G, 


_ 27r (27r)3/2 

= “ (47r/3)H3i?TH 


as the wavenumber correspond- 


( 45 ) 
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The particular value of Rq (or -Rth) appropriate for use in calculating the mass fraction 
which collapses out of the background is uncertain. Our understanding of the detailed 
conditions necessary for the formation of planets and intelligent life has not yet advanced to 
the point of determining what the minimum mass condensation is which is capable of forming 
astronomers. Roughly speaking, we should hlter out condensations that are too small to 
retain metals produced in the hrst generation of stars. The minimum mass condensation 
which is capable of this is currently unknown, however. It is not yet established, for example, 
whether globular clusters of mass 10 ® — lO®Af 0 are capable of self-enrichment, whereby a 
hrst generation of stars generates and releases heavy elements without expelling them from 
the cluster, so that they can subsequently be incorporated in a second generation of stars. 
Dwarf galaxies of even greater mass, in fact, are often postulated to undergo an initial burst 
of massive star formation which leads to supernova-driven expulsion of their interstellar gas 
(containing heavy elements). Even the typical galaxy in a rich cluster of galaxies is widely 
believed to have released most of its heavy elements into the intracluster medium, in order 
to account for the nearly solar metallicity of that gas, which dominates the baryonic mass of 
the cluster. In short, we do not currently know what the minimum mass scale (or associated 
wavelength of density huctuations) is which satishes the necessary condition that the metals 
produced by the hrst generation of stars are retained. Nor do we know if this is a sufficient 
condition for the formation of astronomers. 

In fact, all we can say with certainty is that our own Milky Way galaxy met the necessary 
and sufficient conditions for forming planets, life, and astronomers. The Milky Way has a 
luminosity which makes it roughly an iv*-galaxy, the characteristic luminosity in the bright 
end of the galaxy luminosity function. If the minimum mass scale Mf that can be responsible 
for astronomer formation corresponds to that of an iv*-galaxy, then the data on the galaxy 
luminosity function and the mass-to- light ratio of the bright inner parts of held spiral 
galaxies yields an estimate of the baryonic mass of the bright inner part of an iv*-galaxy of 
Mf ^Mq (see, for instance, Peebles 1993, pp. 122-123). This leads to an estimate of 

i?G — h“®/^(OB/i^/0.015)“®/^Mpc in present units, assuming a cosmic mean baryon density 
which is consistent with the current big bang nucleosynthesis abundance constraints. If, on 
the other hand, we take Afj to be the mass of all the baryons initially within a comoving 
sphere whose volume equals that which, on average, typically contains just one iv*-galaxy 
today, this gives Rq = 2h~^ Mpc. In view of the fact that the Milky Way is actually not 
an isolated galaxy, but is, instead, a member of the Local Group of galaxies, which includes 
more than one iv*-galaxy, we might even wish to consider the possibility that galaxy group 
membership is somehow essential to the formation of astronomers.®'^ In that case, a value as 


^■^For example, group membership might enable a galaxy which undergoes an early burst of star formation 
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large as Rq ~ 3h ^Mpc would even be reasonable. 

In what follows, we will consider a range of values of i?G, therefore. For the case where 
Mf corresponds to the bright inner part of an iv*-galaxy, we shall take Rq = 1 Mpc in present 
units (where, for simplicity, we shall drop the weak dependence of Rq on h for a hxed value 
of We will also bracket the range of possible outcomes by taking values of Rq which 

are smaller and larger than this, respectively. On the low side, we take Rq = 0.01 Mpc, 
relevant for instance, if Afj corresponds to the mass of a globular cluster. On the high side, 
we take Rq = 2 Mpc or 3 Mpc to illustrate the possibilities that Mj corresponds either to 
the total mass within the mean volume per iv*-galaxy or else the mass of a small group of 
galaxies, respectively. 


4.2. The Cold Dark Matter Model 

In order to evaluate a for a given value of i?G, we must adopt a model for the density 
fluctuation power spectrum at recombination. The cosmic microwave background anisotropy 
measured at large angles by the CODE satellite is consistent with Gaussian random noise 
density fluctuations with a scale-invariant primordial power spectrum P{k) oc A;”, where 
n = 1, the case referred to as the Harrison-Zel’dovich spectrum. The range currently allowed 
by a statistical analysis of the hrst four years of data from the CODE DMR experiment is, in 
fact, n = 1.2 ± 0.3 (Bennett et al. 1996). In what follows, we shall generally assume n = 1, 
which is the standard prediction of inflationary cosmology. Later, we can consider the effect 
of a “tilt” in the primordial spectrum away from the shape for n = 1. (Values of n < 1 can 
result, for example, if the primordial fluctuations include a gravitational wave contribution.) 

In general, the power spectrum at recombination differs from the primordial shape A;”, 
except in the long wavelength limit measured directly by COBE. The difference reflects the 
linear growth of the density fluctuations prior to the recombination epoch, which is different 
for different wavelengths. The best-studied and most successful model for the growth of 
density fluctuations to date is the Cold Dark Matter (CDM) model. This model treats 
the CDM density fluctuations as adiabatic fluctuations in a cold, pressure-free gas. Since 
we are interested in the growth of density fluctuations in the baryon-electron fluid, which 
must be present to form stars, planets, and people, we make the assumption that this 
component collapses out in lock-step with the dark matter component, at least for density 


to expel its heavy elements into the surrounding intra-group environment. The latter might then act as a 
reservoir from which the galaxy could later accrete some of its lost metals, after the expelled gas has cooled 
off. 
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fluctuations which are of wavelength large enough to behave in a pressure-free manner. 
As long as we restrict our attention to the epoch of recombination and later epochs and 
to wavelengths larger than the baryon Jeans length in the intergalactic medium, that is, 
the CDM and baryon power spectra should be identical. [For a detailed discussion of the 
effects of Jeans-mass- Altering on the linear growth of baryon density fluctuations in a flat, 
matter-dominated CDM model in which the Jeans mass is increased by the reheating of the 
intergalactic medium which accompanies its reionization, the reader is referred to Shapiro, 
Giroux, & Babul (1994)]. 


4 . 2 . 1 . The Power Spectrum 

We use for the CDM power spectrum the expression given by Fiddle et ah (1996) and 
references therein: 

P{k.z)=2r'‘{4^f*"(S„Yk’'T'‘{q)A-\z.O). (46) 


where 


-4(.-.0) 


'■+(») 
*+(*) ’ 


(47) 


Tiq) = [1 + 3.896? + (W.lqY + (5.466?)^ + (6.716?)"] , (48) 


k 

^ “ hr Mpc-" ’ 

r = , (50) 

and Sff is the dimensionless amplitude at horizon crossing, which must be taken from obser¬ 
vations of anisotropies in the microwave radiation background; Oq is the total matter density 
parameter (Oq = PoIPcritfl^ where Pcrit,o = ISttG); (Ibo is the corresponding parameter 
for baryonic matter; Hq is the Hubble constant; h = iJo/100kms“"Mpc“"; (f_(_ is the pure 
growing mode solution for the evolution of linear density fluctuations in this flat universe 
with nonzero cosmological constant; and A{z,D) = (f_(_( 0 )/(f+( 2 ;) is the linear growth factor 
between redshift 2 ; and the present. For n = 1, these formulae describe the case of the 
Harrison-Zel’dovich scale-invariant power- law primordial spectrum, modified by the growth 
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of fluctuations in a CDM model universe, for a flat universe with a nonzero cosmological 
constant Aq = 1 — Oq = pv jPcvitfi- The htting formula for T{q) is from Bardeen et ah (1986), 
but with r given by a ht by Sugiyama (1995) in the form quoted by Liddle et ah (1996). 
(The numerical coefficients in this formula depend on the present microwave radiation en¬ 
ergy density and on the quantities 100 km/sec and 1 Mpc used in dehning the dimensionless 
quantities q and h, but not on po or Hq.) The formula for T given in equation (50) includes 
an exponential correction factor for the effect of nonzero baryon density. Since the variable 
T is often used in the literature to refer, instead, just to the product ftoh (i.e. without the 
exponential correction factor), the so-called “shape parameter” for CDM models, we will 
also dehne To = Ooh, to use whenever we wish to refer only to this product. In the following 
calculations, we use O^o = 0.015/i“^, consistent with big bang nucleosynthesis constraints 
from the abundance of light elements (e.g. Copi, Schramm, & Turner 1995). 

It is important to note that the explicit dependences of equations (46)-(50) on the 
values of Oq and h do not mean that the power spectrum at recombination is different for 
different subuniverses with different values of the cosmological constant. All factors which 
depend upon Oq and h reflect the fact that a knowledge of the local values of Oq and h in 
our own subuniverse is required in order to interpret present-day observations in our own 
subuniverse (such as those of cosmic microwave background anisotropy) unambiguously to 
determine the power spectrum at recombination assumed common to all subuniverses. We 
must, therefore, distinguish clearly between the particular values of Oq = 1 ~ Aq and h in our 
own subuniverse, on which our inference of the universal power spectrum depends, and the 
variables Oq = 1 ~ Aq and h, different for different subuniverses, on which the probability 
of galaxy formation in any subuniverse depends. To avoid any possible confusion on this 
point, we will, henceforth, indicate the values of these quantities in our own subuniverse by 
adding an asterisk to the symbol (i.e. Aq, Oq, Hq, py, etc.). This notation is not necessary 
for the quantities p, ps, P{k), or cr, however, since these are assumed not to vary from one 
subuniverse to another. 


4-2.2. The Linear Growth Rate in a Flat Universe with Nonzero Cosmological Constant 

In order to evaluate equation (46) for the CDM power spectrum for any particular value 
of the vacuum energy density, we must evaluate the growth factor A(z, 0; Aq) in equation (47) 
as a function of 2 ; and Aq = pv jPent,o- It is convenient to express this in term of a function 
f{\o,z), dehned as the ratio of the growth factor (1 + z) in an Einstein-de Sitter universe 
to A(z, 0; Ao): 
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1 + -2 / N ^+{z) 

/(Ao, z) = = (1 + 2 ;)- 


(51) 


A(z, 0 ) ' '^^^+( 0 )’ 

where is the amplitude of the linear growing mode, which is given for general Aq 7 ^ 0 by 


(Martel 1991), with 


fl \l/2 p fjiu 

S4z-,\„)=(- + l) „,.,6(l + „,)3/2 


Pv ^0 , x_3 

?/ = ^ = 7r(i + ^) 


p{z) Oo 

and Oo = 1 ~ Aq. Using equations (51)-(53), we get, after some algebra. 


(52) 


(53) 


/(Ao, z) — 0,J (1 + z) 


n5/2 


1 + 


dw 


An 


1/2 


nAo/fio 


dw 


Oo(l + zY\ [Jo n;^/®(l + n;)3/2 


-1 


/ no(i + 2 )^ _ 

/o + n;)3/2 


(54) 


For 1 + 2 ; 1, this gives the ^(-independent result 


/(Ao,^)~/(Ao) 


GAq^® ' dw 

50^^ --^0 n;^/®(l -|- 


(55) 


The corrections are of order (l + 2 ()“^, which for the case 2 ( sa 1000 that interests us is entirely 
negligible. We have evaluated the integral in equation (55) numerically, and have plotted 
the function /(Aq) in Figure 4. As we see, /(Aq) differs substantially from unity only for 
relatively large values of Aq. 


4-2.3. Normalization at Recombination from Cosmic Microwave Background Anisotropy 

Measurements 

According to Bunn & White (1996), the hrst four years of data on the cosmic microwave 
background temperature anisotropy detected by the COBF DMR experiment may be ht with 
a dimensionless amplitude at horizon crossing given by the formula 

s*^ = 1.94 X io-5(0*)-o-785-o.o5infi* - 1) + b{n - 1)2] , 


(56) 
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Fig. 4.— Ratio /(Aq) of the linear growth factors for the Einstein-de Sitter model and the 
flat Ao 7 ^ 0 model dehned by equation (55), versus Aq. 

where an asterisk, recall, denotes that the quantities are evaluated for our own subuniverse 
only. There are two sets of values of the constants a and 6, which correspond to the cases 
of n 7 ^ 1 without any gravitational wave contribution (a = —0.95, b = —0.169) and of 
power-law inflation with gravitational waves (a = 1, 6 = 1.97), respectively. 

Equations (46)-(51), along with equations (55) and (56), can now be evaluated to com¬ 
pute the power spectrum at recombination for any flat model for any values of Aq (or, 
equivalently, Oq) and h*. The results for n = 1 are shown in Figure 5 for Aq = 0 and 
h* = 0.5 or 1 (top panel), and for various values of Aq and h* = 0.5 (bottom panel). 


4.3. Results for a and a^p 

The variance cr^ at recombination is given by equations (38) and (46)-(51), as 


h>^+3)/2 


E("+3)/2h^A(zrec,0)-'/Ey2((?„ 


= (ciOo)("+')/"(l + ^rec)-' f {Xo) 


(57) 




log [k (Mpc )] 


Fig. 5.— CDM power spectrum at recombination (^Tpec = 1000), for n = 1, versus comoving 
wavenumber k (= 27r/A, where wavelength A is in present units of Mpc) . Top panel: 
Einstein-de Sitter model (Oq = 1, Aq = 0) with Hubble constants h* = 0.5 and 1. Bottom 
panel: Flat models with h* = 0.5 and A^ = 0, 0.2, 0.4, 0.6, and 0.8. 
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where the symbol labelling the brackets in equation (57) and in what follows, indicates 
that all quantities inside the brackets are evaluated using the values Aq and h* in our own 
subuniverse; Cioo = 2997.9 is the speed of light in units of 100 km/sec; the second equality 
refers to the result to leading order in (1 + and 


(^n 


= < 


/ q'^+‘^T\q)WM2% 

^0 ^ 9max 


9 




dq 

(27r)3/2 q 


L(47r/3)i/3 q^ 


(Gaussian) 
dq (Top-Hat) 


(58) 


The integrals in equation (58) are evaluated numerically. The results are shown in Figure 6 
including the values n = 1, 0.9, and 0.8. (Note: For n 1, we hereafter adopt constants a 
and b for the case of n 1 with no gravitational waves. The case with gravitational waves 
yields a slightly smaller value of a for the same value of n.) 



Fig. 6.— The dimensionless integral AG(<?max) vs. (j'max, dehned by equation (58) for CDM 
density fluctuations, for the Gaussian (solid) and Top-Hat (dashed) window functions, re¬ 
spectively, for n = 1, 0.9, and 0.8, as labelled. 

It is customary to report the normalizations of the power spectra for different models 
in terms of the value of a evaluated at the present for a particular hlter scale, assuming 
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fluctuations continue to grow at the linear growth rate. For our purpose here, however, we 
must evaluate a at recombination, thereby undoing the effects of the growth of fluctuations 
since that epoch which influence the value of a at the present. Our goal, recall, is to use 
the observations of the cosmic microwave background anisotropy made by astronomers in 
our own subuniverse to infer the universal density fluctuation distribution, common to all 
subuniverses at Unfortunately, our ability to infer this universal density fluctuation 

distribution is limited by the fact that we must know the values of Aq and h* in our own 
subuniverse in order to interpret the cosmic microwave background anisotropy measurements 
unambiguously. We illustrate the dependence of the inferred density fluctuations on the as¬ 
sumed values of Aq and h* in Figure 7, where we have plotted the value of a at recombination 
as a function of Aq, for Rq = 0.01,1,2, and 3Mpc, for n = 1, 0.9, and 0.8. We have taken 
Zj-pc = 1000 for the results in Figure 7. The effect of a “tilt” to n < 1 is to decrease a relative 
to its value for n = 1, for the same value of (j'max, or equivalently, of Rq. 

As already mentioned, results for the probability distribution of pv actually depend not 
on cr or p, but on the parameter a^p. The total matter density p at recombination is related 
to the present matter density po by p = po(l + -^rec)^, so as promised cr^p is independent of 
the precise value chosen for Zyf.c'- 

w)l' • (59) 

J * 

5. RESULTS 

With all of the ingredients necessary to evaluate the probability distribution Vohs{pv) 
thus assembled, we can now evaluate the probability of observing any particular value of pv 
anywhere in the universe, as well as the average and median values observed, as functions 
of the values we adopt for py (or Aq = 1 — Oq) and h* in our own subuniverse. There 
are two ways to use this information. We may try to guess the actual value of py, by 
assuming that we live in a typical subuniverse, in which py is equal to either the mean or 
the median observed values of py for all astronomers in all subuniverses. (This will be done 
in subsections 5.1 and 5.2.) Such a “prediction” carries low confidence, since it is always 
possible that py in our subuniverse could be significantly different from the anthropic mean or 
median. Alternatively, we can use what data we have to estimate the range of observationally 
allowed values of py, and then calculate the likelihood that a randomly chosen astronomer 
in any subuniverse would And such a value. (This will be the subject of subsections 5.3 and 
5.4.) If this likelihood proves to be appreciable, then the anthropic principle would survive 
as a possible explanation for the particular value of the cosmological constant in our own 



Fig. 7.— (a) The rms density fluctuation at recombination (i.e. at z^ec = 1000), a = 
(Tree, in the COBE-normalized, flat CDM model, versus Aq, for i?G = 0.01,1,2, and 3Mpc, 
respectively, as labelled, for h* = 0.5, for n = 1 (top panel), 0.9 (middle panel), and 0.8 
(bottom panel), (b) Same as Fig. 7(a), except for h* = 1. 
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subuniverse. 


5.1. The Average Observed Value of the Vacuum Energy Density in a 

COBE-Normalized CDM Universe 

From equations (24) and (59), the average observed vacuum density is given by 


< pv > 
Po 


625(27r)i/2 


243 


Iv 


Io{s) 




(60) 


where again the symbol * indicates that all quantities inside the braces are evaluated at the 
particular values Aq = Aq, Oq = h = h* for our own subuniverse. 

We have used equation (60) to plot the ratio Pvl < pv > versus the assumed value Aq of 
the cosmological constant in our subuniverse, and also versus the ratio Aq/Oq, in Figure 8 
for h* = 0.5 and f, Rq = O.Of, f, 2, and 3 Mpc, and n = 1, 0.9, 0.8. 

For what values of Xq and h* is the vacuum energy density in our own subuniverse equal 
to the observed average of all the subuniverses? For the particular values of Rq and h* 
assumed in plotting the curves in Figure 8 , the intersection of the horizontal dashed line 
with the curves indicates the values for which pf =< pv >■ To answer this question more 
generally as a function of Rq and h*, we have solved the implicit equation pf =< pv > 
numerically by setting the quantity (60) equal to pflPo = ^o/(l ~ ^ 0)5 using the secant 
method. The results are shown in Table 2 for n = 1 and s = f. According to these results, if 
the probability that observers are created in any subuniverse is proportional to the amount 
of mass which eventually collapses out into the bright inner parts of iv*-galaxy-mass objects 
or larger objects (i.e. Rq ~ IMpc), then our own subuniverse has the average observed 
value of the vacuum energy density if the value of Aq which we observe locally is 


Pv 

* 

Po 


M.f , 

1 

r 0.80, 

h* 

= 0.5 ; 


< 7.3 , 

\ * rsj J 
"^0 — 1 

0 . 88 , 

h* 

= 0.7 ; 

(61) 

. 12.4 , 

1 

[ 0 . 93 , 

h* 

= f . 



If the relevant collapsed fraction is, instead, that which condenses into objects as large as or 
larger than those which contain the mean total baryon mass per iv*-galaxy (i.e. Rq ~ 2 Mpc), 
then our own universe has pf =< pv > ii 
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Fig. 8.— (a) The vacuum energy density adopted for our own subuniverse, py = 
‘i{HQ Y\qIG, is plotted in units of the mean value of the vacuum energy density ob¬ 
served in all subuniverses, < py >, versus Aq, the cosmological constant in our own sub¬ 
universe, for local Hubble constant h* = 0.5, for a COBE-normalized flat CDM model with 
primordial power spectrum index n = 1 (top panel), 0.9 (middle panel), and 0.8 (bottom 
panel), assuming shape parameter s = 1. Intersections of the dashed horizontal lines at 
Py! < pv >= 1 with the various curves give the solution of the implicit equation obtained 
by setting < py >= Py in equation (60); (b) Same as Fig. 8(a), except plotted versus Aq/O)); 
(c) Same as Fig. 8(a), except for h* = 1; (d) Same as Fig. 8(b), except for h* = 1. 





















Table 2. LOCAL VACUUM ENERGY DENSITIES WHICH EQUAL THE GLOBAL 

AVERAGE {p*y =< pv >y 


Rg ( Mpc )2 


h * = 0.5 



h * = 0.7 



h * = 1 


* O 

o* 

ILq 

O * 

o * 

* O 

o* 

ILq 

O * 

o * 

* O 

o* 

ILq 

O * 

o * 

0.003 

0.8905 

0.1095 

8.1324 

0.9373 

0.0627 

14.949 

0.9635 

0.0365 

26.397 

0.010 

0.8839 

0.1161 

7.6133 

0.9327 

0.0673 

13.859 

0.9602 

0.0398 

24.126 

0.030 

0.8740 

0.1260 

6.9365 

0.9256 

0.0744 

12.441 

0.9562 

0.0438 

21.831 

0.100 

0.8593 

0.1407 

6.0173 

0.9161 

0.0839 

10.919 

0.9500 

0.0500 

19.000 

0.300 

0.8397 

0.1603 

5.2383 

0.9032 

0.0968 

9.3306 

0.9415 

0.0585 

16.094 

1.000 

0.8046 

0.1954 

4.1177 

0.8791 

0.1209 

7.2713 

0.9252 

0.0748 

12.369 

2.000 

0.7712 

0.2288 

3.3706 

0.8554 

0.1446 

5.9156 

0.9073 

0.0927 

9.7875 

3.000 

0.7418 

0.2582 

2.8730 

0.8339 

0.1661 

5.0205 

0.8912 

0.1088 

8.1912 

6.000 

0.6627 

0.3373 

1.9647 

0.7669 

0.2331 

3.2900 

0.8326 

0.1674 

4.9737 

10.000 

0.5424 

0.4576 

1.1853 

0.6408 

0.3592 

1.7840 

0.6680 

0.3320 

2.0120 

15.000 

0.3582 

0.6418 

0.5581 

0.3539 

0.6461 

0.5478 

0.1665 

0.8335 

0.1998 

20.000 

0.1861 

0.8139 

0.2287 

0.1217 

0.8783 

0.1386 

0.0377 

0.9623 

0.0918 


^For Harrison-Zel’dovich scale-invariant primordial power spectrum (n = 1 ), and s = V/U = 1 . 
^Present units 
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' 3.4, 

1 

ro.77. 

h* = 0.5 ; 


< 5.9, 

\ * rs-( J 

"^0 — i 

0.86, 

h* = 0.7 ; 

(62) 

. 9.8, 

1 

i 0.92, 

h* = l . 



If, instead, galaxy groups are the minimum scale of interest (i.e. Rq ~ 3Mpc), then our 
own subuniverse has py =< pv > H 


' 2.9, 

1 

[0.74, 

h* = 0.5 ; 


< 5.0, 

\ * rs-( J 

"^0 — i 

0.83, 

h* = 0.7 ; 

(63) 

. 8.2, 

1 

i 0.89, 

h* = l . 



On the other hand, if globular cluster formation is enough to satisfy the anthropic constraint 
(i. e., Rg ~ 0.01 Mpc), then our own universe has py =< pv > if 


Po 


' 7.6, 

1 

[ 0.88, 

h* 

= 0.5 ; 


< 13.9, 

\ * rs-( J 

"^0 — i 

0.93, 

h* 

= 0.7 ; 

(64) 

. 24.1, 

1 

i 0.96, 

h* 

= 1 . 



We have assumed s = 1 in all these cases. 


The effect of the “tilt” from n = 1 to n < 1 on these results is to reduce the values 
of < Pv >, or < Ao/Oo >, somewhat. This can be understood in terms of the ratio cr^/cr^, 
where the subscript refers to the index of the primordial power spectrum. This ratio is given 
by 


a 

a 


3 

— 

3 

1 


PcV*h* 


Mpc 


-1 


3(ra-l)/2 



3 

(^max) 



- 1 (^max) - 


3/2 


— (ciooh*)^^” exp [3a(n — 1) + 36(n — 1)^] 






3/2 


(65) 


The exponential factor in the second line of equation (65) is always close to unity for 1 — n ^ 
1. The hrst factor in this line, however, decreases as (1 — n) increases. For (/max ^ 10, the 
results plotted in Figure 6 show that the ratio KnjKi also decreases as 1 — n increases, albeit 
slowly for 1 — n ^ 1. For F* = 0.2, h* = 0.5, and Rq = 1 Mpc (i.e. (/max — 63), for example, 
we hud that cr^jcrf = 0.416 and 0.173 for n = 0.9 and 0.8, respectively. A small tilt to n < 1, 
therefore, decreases the variance of the density fluctuations on galaxy-mass scales. For a tilt 
within the range allowed by the cosmic microwave background anisotropy measured at large 
angles by the COBF satellite, this effect of the tilt translates into a modest decrease in the 
average py observed for all subuniverses. For example, our own subuniverse has the average 
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observed value of pv ( 

he. py =< pv >) 

for Rq = 

: 1 Mpc and h 


1 

r 0.827, 

n = 0.8; 


\ * rsj J 

"^0 — i 

0.857, 

n = 0.9; 



[ 0.879, 

n = 1. 

If i?G = 2 Mpc and h* 

= 0.7, this becomes 



1 

r 0.792, 

n = 0.8; 


\ * rsj J 

"^0 — i 

0.828, 

n = 0.9; 


1 

i 0.855, 

n = 1. 


( 66 ) 


(67) 


Again, we have assumed s = 1 for all these cases. 


5.2. The Median Observed Value of the Vacuum Energy Density in a 

COBE-Normalized CDM Universe 


We can calculate the median value of the vacuum energy density observed in all sub¬ 
universes, (/0v)i/2, by using the solutions for ai /2 = {pv)il2l < Pv > given in equation (32), 
multiplied hy < pv > as calculated already from equation (60). The latter depends upon 
the local values adopted for Aq and h*, used in calculating the variance cr^. Since the value 
of ai /2 is, for a given s, just a number, independent of cr, it is sufficient for us to write 
{pv)i/2 = cti/2{pv) (as for instance ai /2 = 0.43411 for s = 1) and use our previous results 
for < Pv >■ For convenience, we have plotted Pvj{pv)il 2 versus Aq and versus Aq/Oq, for 
h* = 0.5 and 1, Rq = 0.01,1, 2, and 3 Mpc, and n = 1, 0.9, 0.8 in Figure 9. 


For what values of Xq and h* is the vacuum energy density of our own subuniverse equal 
to the median value {pv)i /2 observed in all subuniverses? For the particular values of Rq and 
h* assumed in plotting the curves in Figure 9, the intersection of the horizontal dashed lines 
with the curves indicates the values for which pf = {pv)i/ 2 - To answer this question more 
generally as a function of Rq and h*, we must solve the implicit equation {pv)i /2 = Pv- This 
is similar to our previous implicit equation < pv >= Py, except replaced by the equation 
®i /2 < Pv >= Pv, where ai /2 is a constant for a given s (and only weakly depends on s). 
The results are shown in Table 3 for n = 1 and s = 1. We hud that for Rq = 1 Mpc, our 
own subuniverse has a vacuum energy density equal to that of the median if 


t± 

Po 


' 3.4 , 

1 

[0.77, 

h* 

= 0.5 

< 6.1 , 

\ * rsj J 

"^0 — i 

0.86, 

h* 

= 0.7 

. 10.3 , 

1 

io.91. 

h* 

= 1 ; 


rsj 


( 68 ) 
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Fig. 9.— (a) The vacuum energy density adopted for our own subuniverse, py, is plotted in 
units of the median value of the vacuum energy density observed in all subuniverses, (/ 0 v)i/ 2 , 
versus Aq, the cosmological constant in our own subuniverse, for local Hubble constant 
h* = 0.5, for a COBE-normalized, flat CDM model with primordial spectral index n = 1 
(top panel), 0.9 (middle panel), and 0.8 (bottom panel), assuming shape parameter s = 1. 
Intersection of the dashed horizontal lines at Pv j{pv)il 2 = 1 with the various curves give the 
solution of the implicit equation obtained by setting {pv)il2 = Pv equations (27), (28), 
and (59); (b) Same as Fig. 9(a), except plotted versus Aq/O)); (c) Same as Fig. 9(a), except 
for h* = 1; (d) Same as Fig. 9(b), except for h* = 1. 




















Table 3. LOCAL VACUUM ENERGY DENSITIES WHICH EQUAL THE GLOBAL 

MEDIAN {p*y = {pv)i/2y 


Rg ( Mpc )2 


h * = 0.5 



h * = 0.7 



h * = 1 


* O 

o* 

ILq 

O * 

o * 

* O 

o* 

ILq 

O * 

o * 

* O 

o* 

O * 

o * 

0.003 

0.8795 

0.1205 

7.2988 

0.9301 

0.0699 

13.306 

0.9592 

0.0408 

23.510 

0.010 

0.8709 

0.1291 

6.7459 

0.9239 

0.0761 

12.141 

0.9553 

0.0447 

21.371 

0.030 

0.8592 

0.1408 

6.1023 

0.9164 

0.0836 

10.962 

0.9504 

0.0496 

19.161 

0.100 

0.8414 

0.1586 

5.3052 

0.9049 

0.0951 

9.5152 

0.9429 

0.0571 

16.513 

0.300 

0.8175 

0.1825 

4.4795 

0.8891 

0.1109 

8.0171 

0.9324 

0.0676 

13.793 

1.000 

0.7740 

0.2260 

3.4248 

0.8590 

0.1410 

6.0922 

0.9115 

0.0885 

10.299 

2.000 

0.7319 

0.2681 

2.7300 

0.8273 

0.1727 

4.7904 

0.8884 

0.1116 

7.9606 

3.000 

0.6940 

0.3060 

2.2680 

0.7988 

0.2012 

3.9702 

0.8659 

0.1341 

6.4571 

6.000 

0.5877 

0.4123 

1.4254 

0.7036 

0.2964 

2.3738 

0.7761 

0.2239 

3.4663 

10.000 

0.4279 

0.5721 

0.7480 

0.5080 

0.4920 

1.0325 

0.4584 

0.5416 

0.8464 

15.000 

0.2231 

0.7769 

0.2872 

0.1889 

0.8111 

0.2389 

0.0760 

0.9240 

0.0823 

20.000 

0.0972 

0.9028 

0.1077 

0.0582 

0.9418 

0.0618 

0.0179 

0.9821 

0.0182 


^For Harrison-Zel’dovich scale-invariant primordial spectrum (n = 1), and s = V/U = 1. 
^Present units 
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while if i?G = 2 Mpc, these values shift to 



' 2.7, 

' 0.73, 

h* = 0.5 ; 


Pv ~ ^ 

4.8, A* = < 

0.83, 

h* = 0.7 ; 

(69) 

Po 

. 8.0, 

. 0.89, 

h* = l . 


If i?G = 3 Mpc, then our subuniverse has the median observed value of pv if 


* 

' 2.3, 

' 0.69, 

h* = 0.5 ; 


Pv ~ ^ 

4.0, A* = < 

0.80, 

h* = 0.7 ; 

(70) 

Po 

. 6.5, 

. 0.87, 

h* = 1 . 


If, on the other hand, Rq = 

: 0.01 Mpc, then this happ' 

ens if 


* 

•6.7, 

' 0.87, 

h* = 0.5 ; 


Pv ~ 

12.1, A* = ■ 

0.92, 

h* = 0.7 ; 

(71) 

Po 

, 21.3, 

. 0.96, 

h* = 1 . 


Here again, the effect 

of a tilt to n < 

1 is to 

decrease the values of py 

and of Aq 


corresponding to the median pv somewhat. For example, for Rq = 1 Mpc and h* = 0.7, our 
own subuniverse has the median value of pv (assuming s = 1) for 


For i?G = 2 Mpc and h* 


1 

r 0.80, 

n = 0.8 ; 


\ * r\j \ 

"^0 — i 

0.83, 

II 

o 

(72) 


i 0.86, 

n = 1 . 


0.7, this becomes 



1 

ro.75. 

n = 0.8 ; 


\ * rsj J 

"^0 — i 

0.79, 

II 

o 

(73) 

1 

i 0.83, 

n = 1 . 



5.3. Observational Constraints on the Cosmological Constant in Our Own 

Subuniverse 

Ostriker & Steinhardt (1995) have argued that the apparent discrepancies that had 
previously been identihed between observations and the predictions of the standard CDM 
model (i.e. cold dark matter in a flat universe with zero cosmological constant and the scale- 
invariant Harrison-Zel’dovich primordial power spectrum) can be reconciled if the standard 
CDM model is modihed to admit a nonzero cosmological constant roughly in the range 
Aq = 0.65 it 0.1. We have reproduced the observational and theoretical constraints that led 
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h* 

Fig. 10.— (a) Observational Constraints on Aq and h* for own subuniverse. Curves labelled 
“LSS” and “Fq = 0 .2” or “Fq = 0.3” bound the region allowed by the constraint Fq = Oq^* = 
0.25 it 0.05 derived by matching the spatial and angular correlation statistics from galaxy 
surveys with the theoretical predictions of the large-scale clustering of galaxies in a COBF- 
normalized, flat CDM model with primordial power spectrum index n = 1. The curves 
labelled “(Js X-ray clusters” bound the values of Aq and h* which make this CDM model 
satisfy the constraint on the present space density of X-ray clusters. The curve labelled 
“to = 12 Gyr” indicates the lower limit which makes the age of the universe at least as large 
as current estimates of the minimum age of globular clusters. The curves labelled 
indicate the boundaries dehned by the X-ray-measured total and baryonic masses of clusters 
of galaxies, together with the big bang nucleosynthesis limits on the baryon mean density 
and the assumption that the ratio of baryon to total mass inside each cluster equals the ratio 
of universal mean values. The curve labelled “gravitational lenses” indicates the upper limit 
imposed by counts of quasars lensed by intervening galaxies. The dashed curves labelled 
“(py)i/ 2 ” are the values for which our own subuniverse has the median value of pv for all 
subuniverses, if Rq = 1 Mpc, if n = 1 (top dashed curves), 0.9 (middle dashed curve), or 0.8 
(bottom dashed curve), (b) Same as Fig. 10(a), except Rq = 2 Mpc. (c) Same as Fig. 10(a), 
except i?G = 3 Mpc. 
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Ostriker & Steinhardt (1995) to this conclusion, plotted here in Figure 10 as a series of upper 
and lower bounds in the (Aq, /i*)-plane. The constraint curves plotted in Figure 10 are based 
on Ostriker & Steinhardt (1995) and references therein, with the following additions. 

The data on the large-scale clustering of galaxies from galaxy surveys constrain the 
flat, CDM model by requiring that the spatial and angular correlation statistics of the 
observed galaxies in our local universe at the present epoch agree with the predictions of 
structure formation by gravitational instability in the CDM model. This leads to upper 
and lower bounds on the so-called “shape parameter” Fq = Ooh, given by Fq = 0.25 ± 0.05 
(assuming n = 1) (see curves labelled “Fq” and “LSS” in Fig. 10). A similar constraint 
results from the requirement that the CDM model reproduce the observed space density 
and luminosity function of X-ray clusters in the present universe. We plot this constraint 
separately from that of the shape parameter Fq which Ostriker & Steinhardt (1995) plotted. 
This X-ray cluster abundance constraint is expressed by Viana & Fiddle (1996) as bounds 
on the rms density fluctuation (t’^{z = 0) for a smoothing radius (in present units) given by 
i?TH 8h ^Mpc. Assuming n = 1, these bounds are given by 

-(0.59-0.016fi*-|-0.06fi*2)j +32</-(fi0% 

J -244>{U*)% ’ 

where 

(^(0*) = (0*)°-26i°siofio . (75) 

This amounts to a constraint on Aq and h* which is similar to that from the correlation 
statistics. (The above-mentioned bounds from the statistics of large-scale structure (“LSS”) 
in the galaxy distribution also refer to the part of the CDM power spectrum at wavelengths 
A ^ 8/i“^Mpc.) 

Fstimates of the total masses and baryonic mass fractions of individual clusters of galax¬ 
ies, derived by htting the X-ray surface brightness prohles of each cluster and assuming the 
cluster intergalactic medium is an isothermal sphere in hydrostatic equilibrium with a virial- 
ized cluster gravitational potential, yield another pair of bounds. If the assumption is further 
made that the ratio of baryonic mass to total mass of each cluster is equal to the universal 
mean ratio, Qbo/^o, then a comparison of this X-ray-estimated ratio with the constraints 
on Qbo from standard big bang nucleosynthesis and the observed light element abundances 
(i.e. 0.01 ^ flBoi^*)^ ^ 0.02; Copi, Schramm, & Turner 1995) implies a constraint on the 
total density parameter given by 0.09 ^ ^ 0.33. The curves labelled and 

“X-ray cluster masses + big bang nucleosynthesis” in Figure 10 indicate the bounds on Aq 
and h* which result from this argument. Some recent numerical gas dynamical simulations 
of cluster formation in the flat, matter-dominated CDM model suggest that the upper bound 



(74) 
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on which results from the high values estimated for cluster baryonic mass fraction 

by the equilibrium model described above may be too low (e.g. Bartelmann & Steinmetz 
1996; Martel, Shapiro, & Valinia 1996; Valinia 1996). The simulated clusters, when properly 
resolved, are often found to be comprised of subclusters in the act of merging and, together 
with projection effects, this can cause an observer who uses the assumption of isothermal 
spheres in hydrostatic equilibrium to underestimate the total mass and overestimate the 
baryon fraction. 

The estimated minimum age of globular clusters derived by comparison of theoreti¬ 
cal models of stellar evolution with observed globular cluster H-R diagrams is about 12 X 
10®years. This leads to a lower limit to Aq for each h* based on the requirement that the 
age of our universe exceeds this estimate of the minimum age of globular cluster stars (see 
curve in Fig. 10 labelled “to = 12Gyr”). 

An upper bound to Aq results from the comparison of the statistics of quasars observed 
to be gravitationally lensed by intervening galaxies with the predictions of flat cosmological 
models with a nonzero cosmological constant. A flat cosmology with cosmological constant 
tends to produce more gravitationally lensed quasars than does such a cosmology with zero 
cosmological constant. The resulting limit was quoted by Ostriker & Steinhardt as Aq < 0.75. 
More recently, Kochanek (1996) has argued for a somewhat tighter limit. However, limited 
observational data and the possibility that evolution effects on the population of lensing 
galaxies have not been properly taken into account suggest that this limit is still uncertain. 
We plot the upper bound quoted by Ostriker & Steinhardt (1995) in Figure 10, but caution 
that this limit is still in flux and that the future change can be up or down. 

Finally, an independent constraint was recently derived by Perlmutter et ah (1996), 
using observations of Type la supernovae to infer the relationship between redshift and 
distance. Their result is consistent with the limit Aq < 0.5. However, this interesting 
approach is too preliminary to be considered reliable at this time. We have not included it 
in Figure 10. 


5.4. How Likely is the Value of the Cosmological Constant Observed in Our 

Subuniverse? 

Along with the observational constraints discussed in the previous subsection, we have 
plotted in Figure 10 the values of the cosmological constant Aq such that our own universe 
has the median value observed for all subuniverses for different values of h* and for n = 1, 
0.9, and 0.8, for Rq = 1, 2, and 3Mpc and s = 1. It is apparent from Figure 10 that if 



49 


the universe is flat and the inflationary CDM model applies, then the range of values of the 
cosmological constant allowed for our own subuniverse is somewhat below the median, but 
not far from it. 

This statement is sensitive to the choice of the smoothing scale Rq. As discussed in 
Section 4.1, this scale is intended to reflect the minimum size of fluctuations which are 
responsible for forming observers in any subuniverse. The value Rq = 1 Mpc was chosen 
to correspond to fluctuations which encompass roughly the mass of the bright inner part of 
an iv*-galaxy. If, instead, observers are formed if even a much smaller, globular-cluster-size 
object collapses out of the background, then the global median value of pv would be larger 
since the variance of the density fluctuations at recombination Altered on this smaller scale is 
larger. In that case, the median {pv)i /2 is on the high side, somewhat further from the range 
in the (Aq,/ i*)-plane allowed by observations of our own universe. Conversely, if observers 
are formed only if fluctuations larger than 1 Mpc collapse out, then the global median {pv)i /2 
will be smaller [since cr[zj-Qc) is smaller if Rq > 1 Mpc than if Rq = 1 Mpc] and the match 
between the allowed range in the (Aq,/ i*)-plane and {pv)i /2 improves. We illustrate this 
by plotting the median value curves for Rq = 2 and 3 Mpc along with the observational 
constraints, in Figures 10(b) and (c). 

These results also depend weakly on the value adopted for O^q. A small increase of 
O^Q has the effect of decreasing T* and, with it, the value of (t[zj-qc)- This, too, would 
drive the global median {pv)i /2 downward for given values of Rq and h*. We chose O^q = 
0.015(/i*)“^ here, but the big bang nucleosynthesis limits allow a value = 0.02(h*)~^ 
and uncertainties exist even in that upper limit. 

Even though pv in our subuniverse seems to be below the median, the allowed range in 
the (Aq,^*)- plane in Figure 10 includes values that are reasonably likely. To see this, it is 
useful to recast the curves of Figure 10 as curves in the (a*, /i*)-plane, where a = pv j <pv>- 
For each point (Aq,^*) on a curve in Figure 10 we can use the results of Sections 3 and 4 
to compute <py>, and hence the ratio a* = pyj < pv >■ We plot these constraint curves 
in the (a*,/i*)-plane in Figure 11. Inspection of this figure shows that the observationally 
allowed values of a* are between 0.01 and 0.04 for Rq = 1 Mpc, between 0.015 and 0.1 for 
i?G = 2 Mpc, and between 0.03 and 0.15 for Rq = 3 Mpc. For all of these values of i?G, 
there is an appreciable overlap between these observationally allowed values, and the range 
from 0.02 to 3.7 which was found (for s = 1) in equations (33) and (34) to be anthropically 
likely. 
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Fig. 11.— (a) Observational Constraints of Fig. 10 are plotted instead as curves in the 
(a*,/i*)-plane, where a* = pyj < pv >, and < py > is evaluated by assuming COBF- 
normalized density fluctuations for the flat CDM model with n = 1, as inferred by adopting 
the value py = Py, assuming Rq = 1 Mpc and s = 1. Curves are labelled just as in Figure 
10. (b) Same as Fig. 11(a), except Rq = 2 Mpc. (c) Same as Fig. 11(a), except Rq = 3 Mpc. 
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6. SUMMARY AND CONCLUSION 

The range of values of the cosmological constant that allow life to arise is so narrow, 
that within that range, we can assume that the a priori probability distribution of values of 
the cosmological constant is constant. The probability that a particular value of the cosmo¬ 
logical constant is observed in our universe is then proportional to the number of observers 
who might measure that value. That abundance is, in turn, proportional to the fraction of 
matter which eventually collapses out of the background into gravitationally bound concen¬ 
trations large enough to initiate star formation and retain heavy elements, presumed to be 
prerequisites for the origin of planets and intelligent life. We have derived an analytical esti¬ 
mate for this “collapsed fraction” based upon a simple, pressure-free, spherically symmetric, 
nonlinear model for the growth of density fluctuations in a flat universe with arbitrary value 
of the vacuum energy density py, applied in a statistical way to a distribution of cosmolog¬ 
ical density fluctuations. We have evaluated the resulting probability distribution for the 
observed values of pv for density fluctuations which are Gaussian random and of linear am¬ 
plitude at recombination. We hud that the probability distribution in that case is a unique 
function of py/ < py >, where < py > is the average of the observed values of py over 
all subuniverses, and a shape parameter s which characterizes the variation of density with 
radius within a given density fluctuation. The dependence on the parameter s is, moreover, 
found to be very weak. The values of py (such as the median value) at which the integrated 
probability distribution takes any dehnite values are simply proportional to the mean value 
< py > for a given s, with proportionality constants that are fairly insensitive to the value 
of s, and for a given value of s, these values of py are, in turn, exactly proportional to the 
quantity cr^p evaluated at recombination, where is the variance of the density fluctuations 
and p is the cosmic mean matter density. The dependence of this proportionality constant 
on the value of s is, once again, quite weak. The quantity cr^p is assumed to be common to 
all subuniverses in which life can arise, because in these subuniverses py is negligible at and 
before recombination. 

Presumably it will some day be possible to calculate cr^p unambiguously by using as¬ 
tronomical observations to measure the density fluctuations in our own universe. While 
this is not yet possible with great certainty, current measurements of anisotropy at large 
angles in the cosmic microwave background by the CODE DMR experiment represent sub¬ 
stantial progress toward this goal. In particular, the detected cosmic microwave background 
anisotropy hxes the fluctuation amplitude at long wavelengths and is consistent with Gaus¬ 
sian random density fluctuations in a flat universe, with the scale- invariant primordial power 
spectrum, P{k) oc with n sa f, as expected from inflationary cosmology. Unfortunately, 
the precise interpretation of these anisotropy measurements is, itself, dependent upon our 
uncertain knowledge of the values of Aq and h for our observed universe. In addition, the 
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variance cr^ which is relevant to our anthropic probability calculation is that for the density 
fluctuations after the density has been smoothed over some length scale Rq so as to eliminate 
fluctuations which are too small to contribute to the formation of astronomers, and Rq is 
orders of magnitudes smaller than the long wavelengths at which the fluctuations are con¬ 
strained directly by the COBE DMR anisotropy measurements. Unfortunately, therefore, 
neither the value of Rq nor the amplitude of the density fluctuations at wavelengths close to 
i?G are precisely determined at this time. In the meantime, in order to specify cr^p, we have 
here adopted the CDM model for the growth of a scale-invariant primordial spectrum of 
Gaussian random density fluctuations in a flat universe with nonzero cosmological constant, 
with an amplitude set by the COBE DMR data, and parameterized our results in terms of 
the unknown hlter scale Rq and power spectrum index n. 

Our results are encouraging from the point of view of the anthropic hypothesis. Although 
the range of observationally favored values of pv is somewhat less than the median value of 
pv for all subuniverses, it has a signihcant overlap with the range of values between the 5th 
and 95th percentile for all subuniverses. 

In short, these results explain why there might be a nonzero, but small value of pv in 
our universe, if all values of pv are otherwise equally likely to occur. They show that a range 
of values close to those favored by current observations are reasonably probable, while values 
which are orders of magnitude smaller or larger are extremely improbable for us to observe. 
This is the essential ingredient in the anthropic explanation for the observed value of Aq in 
our own universe. 

With a continued improvement in measurements of pv and the spectrum of fluctuations 
at recombination, it may turn out that the actual value of pv in our subuniverse is more 
than two orders of magnitude less than the average. In this case we would have to conclude 
that the anthropic arguments used here do not explain the smallness of the cosmological 
constant. Unfortunately, the converse is not possible; observation of a value of pv that is 
anthropically likely would support the idea that there is a diversity of possible pv values, 
but since we only observe one subuniverse, astronomical observation alone cannot conhrm 
this idea. Ultimately this issue will have to be settled by advances in fundamental physics, 
which we hope will tell us whether in fact it is correct that there are many subuniverses 
with different values of the cosmological constant. If this is not correct, then there is no 
justihcation for the anthropic reasoning used here, while if it is correct, then these anthropic 
arguments are just common sense. 
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